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\ The d-dimensional n-spin facilitated kinetic Ising model is studied analytically starting from 

fSJ . usual master equations and their transformation into a Fock-space representation. The evolution of 

» I ' relevant operators is rewritten in terms of a projection formalism. The obtained frequency matrices 

CIh, and memory terms are analyzed. Especially, the influences of the memory terms is approached by 

■ using standard techniques of the usual mode coupling approach. The temperature dependence of the 

relaxation times related to the n-spin facilitated kinetic Ising model shows a weak non-Arrhenius 
behavior. Furthermore, a characteristic stretched decay of the correlation function is obtained. 
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I. INTRODUCTION 



, In spite of advances in the description of Uquids near the glass transition using different approaches the phe- 
nomcnon is generally not completely understood. Supercooled fluids reveal often a non-Arrhenius behavior of relax- 
^ , ation times and a characteristic stretched exponential decay of correlation functions. But a long range order is not 
^ ■ developed in contrast to conventional phase transitions. The dynamics of the glass transition can be described by an 
^ ' increasing coopcrativity of local processes with decreasing temperature Q . The cooperativity leads to the well known 
' ■ slowing down of the relaxation times which can be illustrated by a strongly curved trajectory in the Arrhenius plot 
(logarithm of the relaxation time r versus the inverse temperature T~^). One possible fit of this curve is given by the 
Q Williams-Landel- Ferry (WLF) relation i.e. Inr c>c (T — To) with a finite Vogel temperature Tq. 
O \ Mode coupling theories (MCT) predict a completely ergodic decay of correlation functions above a critical 

temperature Tc and an incomplete decay below Tc- The incomplete decay is usually identified as the relatively fast 
I /3-process caused by processes related to the breaking up of the local cages. The long time regime of the ergodic decay 
^ , above Tc is denoted as a-process. The relaxation time of the a-process increases rapidly with decreasing temperature 
' and below only the /^-process remains effective, i.e. at Tc the system undergoes a sharp phase transition to a state 
with partially frozen (density) fluctuations. Note that Tc is in the range between the melting temperature Tm and 
' the glass temperature Tg, e.g. > Tc > Tg. 

It is really the a-process exists also below Tc- This process leads to a very slow decay of the apparently frozen 
structures, i.e. the nonergodic structures obtained from the MCT [^J] are approximately stable only for a finite 
[ time interval. This slow decay shows the typical properties which correspond usually to the dynamics of the main 
, glass transition (WLF like behavior of the relaxation time, stretched exponential decay of the correlation function). 
This effects can be partially described in terms of an extended mode coupling theory introducing additional 

^ J hopping processes. There exists also various alternative descriptions which explain the cooperative motion of 

I ' the particles inside a supercooled liquid below Tc- One of these possibilities is the n-spin facilitated kinetic Ising 
■ O ' model [pd|-p^, originally introduced by Fredrickson and Andersen. The base idea of this model consists in a coarse 
graining of space and time scales and simultaneously a reduction of the degrees of freedom. In detail that means: 



1. Coarse graining of spatial scales: The supercooled liquid is separated into cells, so that each cell contains a 
sufficiently large number of particles which realize a representative number of molecular motions, i.e. the many 
body system is considered of a virtual lattice with the unit size I. This lattice has no influence on the underlying 
dynamics of the supercooled liquid. 

2. Reduction of the degrees of freedom: Each cell will be characterized by only one trivial degree of freedom, i.e. the 
cell structure enables us to attach to each cell an observable tTj which characterizes the actual state of particles 
inside the cell j. The usual realization is given by the local density pj (particles per cell) with ctj = — 1 if pj > p 
and Uj = \ ii pj < p where p is the averaged density of the system. This mapping implies consequently different 
mobilities of the particles inside such a cell, i.e. aj = —1 corresponds to a more immobile (or solid like) state 
and aj = 1 to a more mobile (or liquid like) state of the cell j. The set of all observables a ={crj} forms a 
configuration. The evolution of the statistical probability distribution function P (ct, t) can be described by a 
generalized master equation [[l5[ using a projection of the real dynamics onto the dynamics of a: 
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dP (g, t) 
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t-t')P{a, t')dt' 



3. Coarse graining of the time scale: The last step bases on the assumption that possible memory terms K[a, cr', t) 
of the generalized master equation is mainly determined by fast molecular processes while the slow dynamics 
is mainly reflected by the temporally local contributions ^ £'((7, (t')P (cr, i). Of course, the validity of this 

a' 

assumption depends strongly on the choice of the remaining degrees of freedom, and in many cases it is very hard 
(or impossible from the actual point of view) to give a satisfactory explantation of this assumption. However, 
if this separation of the dynamics is justified, an elementary time scale larger than the time scale of the fast 
molecular processes can be introduced. Therefore, the memory will be reduced to simple temporally local terms, 
i.e. if (cr, cr', t — t') = 5 {t — t') K{a, a' , T)dT. One obtains an evolution equation which is equivalent to the 
mathematical representation of a usual master equation. 

^^ = ^L(a,.')P(.,t) (1) 



The dynamical matrix L{a,a') = L'{a,a') + K{(t,(t' ,T)dT is determined by the above discussed formal 
procedure. Unfortunately, a direct calculation is mostly very complicated, so that one should use reasonable 
assumptions about the mathematical structure of L. 

To make the time evolution of t he g lass configurations more transparent we use the argumentation following the 
idea of Fredrickson and Andersen [jllp^, i.e. we suppose that the basic dynamics is a simple (Glauber) process 
aj = +1 <-> CTj = —1 controlled by the thermodynamical Gibb's measure and by self induced topological restrictions. 
In particular, an elementary flip at a given cell is allowed only if the number of the nearest neighbored mobile cells 
{aj = +1) is equal or larger than a restriction number n with < n < Zc (zc- coordination number of the lattice). 
Elementary flip processes and geometrical restrictions lead to the cooperative rearrangement of the underlying system 
and therefore to mesoscopical models describing a supercooled liquid below Tc. These models [^ij-Q are denoted 
as n-spin facilitated Ising model on a d-dimensional lattice SFM[n,d]. The self-adapting environments influence in 
particular the long time behavior of the spin-spin correlation functions and therefore of the corresponding density- 
density correlation functions. The SFM[n, d] was studied numerically [p^-p^ (SFM[2, 2]) and recently also analytically 
@(SFM[1,1]). 

From this point of view it will be an interesting task to derive a set of equations related to the SFM[n, d] which are 
similar to the well known Mori-Zwanzig equations [pl[| and which can be used as a reasonable basis for a further 
treatment, e.g. for a continuous fraction analysis or for a mode coupling approach. 

The aim of the present paper is the derivation of such evolution equations and their analysis in terms of a mode 
coupling approach. We restrict our investigation to the analysis of the SFM[2, d] but a generalization to another class 
of spin facilitated kinetic Ising models is always possible. The starting point is the mapping of the master equations 
of the SFM[2, d] to evolution equations in a Fock-space representation. Using a projection formalism one obtains 
evolution equations for a set of relevant observables and consequently for the corresponding correlation functions. 
The paper ends in an analysis of these correlation functions in terms of the frequency matrices and memory terms. 



II. FOCK-SPACE APPROACH 

Following Doi j2^, compare also |Q, the probability distribution P{a,t) can be related to a state vector \F(t)) 
in a Fock-space according to P{a;t) — {a\F{t)) and \F{t)) — J2a^('^^^) I''')' respectively, with the base vectors |cr). 
Using this representation, the Master equation (^ can be transformed to an equivalent equation in a Fock-space 

d,\Fit))=L\Fit)) (2) 

The dynamical matrix L(cr, a') of (|l|) is now mapped onto the operator L given in a second quantized form with d 
and being the annihilation and creation operators, respectively, for flips processes. Usually L is expressed in terms 
of creation and annihilation operators which satisfy Bose commutation rules |^,^,^ . 

The SFM[n,(i] can be interpreted as a lattice gas (cr^ = 0: empty cell, cr^ — 1: occupied cell) considering the excluded 
volume effect, i.e. changes of the configuration a are possible only under the presence of the exclusion principle. To 
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preserve the restriction of the occupation number in the underlying dynamical equations too, the commutation rules 
of the operators d and d'^ are chosen as those of Pauli-operators ||23|,p3,p6|,p7| : 



[d.,d] = 4,(l-2dld.) [d,,d,] = [dl,4]^0 d^^idl) 



U2 







(3) 



It should be remarked that the method can be extended to the case of higher restricted occupation numbers ||2^ . 
As it was shown by Doi p2[ the average of a physical quantity B{a) is given by the average of the corresponding 
operator B{t) = \a) B{a) ((t| via 



= B F{t) 



(4) 



using the reference state (s| = (a|. The normalization condition is manifested in the relation < s\F{t) >= 1. In 
the same way, arbitrary correlation functions can be expressed by 



A{t)Bit')) - ^ A(a)P(a,i;a',t')S(a') 



Aexp 



{Lit-t')}B\Fit') 



From this point of view, it is possible to create the evolution equations of various averages and correlation functions. 
For example, using (||) and (^) one obtains the evolution equation for an arbitrary operator B j2^ : 



dt(B) = 



[B,L] F{t) 



(5) 



Here we have used the necessary relation < s\L = 0, which is an immediately consequence of the normalization 
condition. The evolution operator for the SFM[2, d] can be written as 



L = Y1 ^mDjDk [f3id, - bi) + X{d\ - (1 - A)) 



(6) 



with the particle number operator Di = d-di and temperature dependent jumping rates A and /?. At^ufc is a lattice 



function with k. 



i\jk 



1 if j fc and j and k are neighbored to lattice cell i. Applying a simple activation dynamics 



one obtains for the jumping rates: 



and 



A i/-i(T)exp(-e/T) 



(7) 



where v ^(T) is an elementary temperature dependent time scale (e is the energy difference between the solid and 
liquid like state). Note that the stationary state corresponds to an average a^q = ^-^j) ~ ^/(^ + /^) which can be 

obtained directly from (||). The knowledge of L and the corresponding evolution equation ^ allows a reasonable 
analysis of the SFM[2, d]. 

The Fock space formalism has the decisive advantage of a simple construction principle for each evolution operator L 
on the basis of creation and annihilation operators. Therefore, this method allows investigations of master equations 
for various evolution processes, e.g. aggregation, chemical reactions [^,^, nonlinear diffusion and just the spin 
facilitated kinetic Ising model. 



III. PROJECTION EQUATIONS 



A. relevant operators 



The dynamics of an arbitrary physical system can be described by a reasonable set of relevant operators. We use the 
normalized local deviations of the configuration from the thcrmodynamical average and the corresponding derivatives 
with respect to the time as suitable relevant observables for the investigation of the SFM[2, d] 



bi{t) - O-pq 



and 



(8) 
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(/3 = 0, 1, ...,5max; usually the upper borderline is a finite integer number, but (/max ^ oo is also possible). These 
covariant operators must be extended by the corresponding contravariant operators 

?7f)(i)=%W and ?,i''\t)=L^m (9) 

Using (H) and (H) we construct the backward projection operator P: 

p-Y.{ v^')9^f^' with E ('^i'^^^^^S' = (10) 

(with a, f3, ...e [0,5max])- The projection operator leads to an identical mapping of the relevant operators onto itself, 
i.e. fff^'^ P — f}"^^ ■ Consequently, the orthogonal projection operator Q is given hy Q = \ — P with "ffjJ^Q = 0. 



B. Basis equations 



The evolution equation (|2|) leads to the formal solution \F{t)) — exp Wi^))- The dependence of \F{t)) on the 
time can be transferred to an arbitrary operator analogous to the transformation of Schrodinger's representation into 
the Heissenberg picture. Therefore one obtains time-dependent operators 



= i?exp|ii| (11) 

re evolution equations for t 

It follows from (hll) 



The derivation of the evolution equations for the relevant observables starts from the formal time evolution of rfjj\t). 



dt 



vV'm (12) 



This equation is the basis for the derivation of projection equations for the relevant observables in analogy to the 
well known Mori-Zwanzig pl| ] equations for classical or quantum mechanical equations of motion. The application of 
1 — P + Q onto the operator L leads to a formal splitting into a relevant and an irrelevant part (Note that P realizes 
a projection onto the subspace L|| of relevant operators, whereas Q projects onto the linearly independent subspace 
L± of all other operators). Hence, 



= 4^)exp{Lt}L = 4^)L(P + g)exp{Lt} 
E^^if +4"^iQexp{Lt} 



(13) 



/3j 

with the frequency matrix 

<' = -T.{^'^'^t')9:f (14) 

The second term of ( [T^ ) can be rewritten by using an identical transformation of exp |^^| iuto an integral expression: 
exp|Li|= / dt' exp ^^LQ{t-t')^ LP exp ^^Lt'^ + exp |LQi| (15) 





This relation allows the derivation of rigorous projection equations similar to the usual Mori-Zwanzig-equations [ |21| , [32[ |: 



E^if + / ^^'E<^(*-o^f (O + /i^^(^) 

n , Jo „ 



0:3 13,3 



(16) 



with the residual forces 
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= V^^\i)LQe^j> [LQt] = /^^) exp [hQt] (17) 
(with the properties f'^\t)Q — fk'\t) and f\?\t) P — 0) and the memory matrix: 

^S'^ it - = E (^i'^^^exp {LQ{t t')} L^i"^) (18) 

The comparison between (|l6|) and the usual Mori-Zwanzig-equations shows a formal equivalence. Both equations con- 
tains frequency terms, memory and residual forces with a similar mathematical structure. But there is a fundamental 
difference which can be studied directly by inspecting the memory kernel. The memory of Mori-Zwanzig-equations can 
be written always as a correlation function of the residual forces. This relation can be interpreted as a representation 
of the fluctuation-dissipation theorem, and it is causally connected with the fact, that the Mori-Zwanzig equations 
are related to reversible classical or quantum mechanical equations. On the other hand, the memory ( p^ cannot be 
completely constructed from residual forces (pT|). The cause is the irreversible character of the underlying master 
equation. 



C. Projection equations for a reduced set of relevant observables 



We choose (/max = 1 for the following investigations, i.e. the relevant observables are fif^\t) = f]i{t) and 
fli^ {t) = fliit)L. This settling corresponds slightly to mechanical systems, which are completely determined by 
spatial coordinates and velocities. The general system of equations (nO) becomes 



i /3=0,1 

9 



-n 



j /3=0,1 L J 



(19) 



(0/3) 



SnjS^I^, fi°\t) = and K^^f\t~t') ^ 0. Therefore, the first 



.(0/3) 



A simple analysis leads to the simplifications 17, 

equation will be reduced to the identity dtf]n\t) — dtf]n{t) — fjn\t) and the second equation can be rewritten as 



9*2 



Vnit) 



E <°)%(t)-fl7:-^^,7,(i) 



(11) 9. 



j 



(20) 



The result is a second order differential equation which reflects the complete dynamics of the relevant observables. 



D. Projection equations for correlation functions 

An important rule for experimental and theoretical investigations of the glass transition plays the time-dependent 
equilibrium correlation functions of the relevant observables. Especially the SFM[n,(i] should be characterized by 

*™n(i) = {fln{t)flrniO)) = |?7„ CXp {ii} t}™ | F(0)^ (21) 

These correlation function is equivalent to the normalized spin-spin correlation: 



0'oq(l — CTcq) 

which should be similar to the normalized density-density correlation of the underlying supercooled liquid, i.e. 
^riTnjt) ^ {6p{r,t)Sp{r' ,0)), see also the above discussed mapping p ^ a. The evolution equation of ^nm{t) follows 
from ( |20| ) by a right hand multiplication with and a subsequent determination of the average. The contributions 

of the residual forces fn'' (t) vanish identically. One obtains a homogeneous integro-differential equation: 
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9^$ 



"J 



9t 



(22) 



This evolution equation of $nm(i) is a rigorous second order integro-difFerential equation, which will be analyzed 
now. To this aim it is necessary to determine the frequency and memory parts for the thermodynamical equilibrium. 
Here, correlation functions, memory and frequency matrices should be homogeneous and isotropic functions, e.g. 
^mn(t) = <i>(|n — m| , <) or Kmn\t) = K'^"'^\\n — m| , i). The homogenity is a direct consequence of the underlying 
translation invariance. On the other hand, isotropy can be expected for the asymptotic case of the continuous limit, i.e. 
for |n — m| oo. However, the isotropy is partially disturbed at finite distances as a consequence of the underlying 
lattice structure. The Fourier transformation of (E3) can be obtained by using the representation 



<i>mn{t) = ^^{(l,t) exp{iq{n-m)} $(q, t) = ^ ^ $(|n|, t) exp {-iqn} 

q 



of the correlation function whereas the frequency matrix and the memory can be written as: 



(23) 



and 



(24) 



respectively, n and m denote the lattice vector of size n and m, q is a vector of the first Brillouin zone corresponding 
to the lattice. We use a cubic lattice for the following calculations but an application of another lattice type is always 
possible. The Mori-Zwanzig-equation ([22|) becomes now: 



^q,t) = -N[n('^Hq)^q,t) + nif{q)<i>iq,t) 

+ N dt' (q, t - t') $(g, t') + if (") {q, t-t')^q,t') 



(25) 



Note that all quantities depend only on q = |q| (at least for the continuous limit) because of the isotropy. Finally, 
the Laplace transformation 



POO 

$(g, z)= / dt exjp{-zt}^{q,t) 
Jo 



leads to the algebraic equation 



4>(g,z) = 



*o(9) 



iVl](l°)(g) - NK(^°\q, z) ^ zgo{q) 
z + Nn^''\q)-NK^''\q,z)+go{q) 



(26) 



(27) 



which considers the initial conditions $(g,0) = $o(<z) and $((7,0) ~ $o('?)- Furthermore, the quantity goiq) denotes 
the ratio goiq) = ^o{q)/^o{q)- It should remarked that especially ri(^^^((?) ^ and K'^^°^q,z) ^ are consequences 
of the irreversible master equations. On the other hand, the usual Mori-Zwanzig equations are founded on reversible 
Liouville operators which lead immediately to r2'^^''(g) = and K'^^^^{q^ z) = 0. 



IV. DETERMINATION OF FREQUENCY MATRICES 

The concrete determination of the frequency matrices is possible by using the concrete evolution operator L (|^). 
Note that the projection equations ( p^ ) are valid for an arbitrary physical system which can be described by master 
equations. 

The frequency matrices contain always the matrix g, which can be calculated from (p^, i.e. g is determined by the 
following system of linear equations: 
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/^r^ = (^r'C9 with y: = (28) 

7, A; 

Using the definition 

rffc = (v^L^f)k) (29) 



one obtains simple expressions for the matrix h (i.e. h}'^^'^^ — r"^^^) as well as for the frequency matrix (p^: 

f^if --Err^'^r (30) 

The knowledge of T^j. {(3 = 0...3) allows the exact determination of Q.^^'^\q) and ^l^^^\q). The values T^^, can be 
obtained straightforwardly by using (^ and the commutation relations (||). It follows: 

= (^-^^ (A"^™„ + S"e„„, + C"Xmn + ^"C«m) (31) 

The lattice functions 0„„i, Xmn and Cnm vanish, except for the following cases: Qnm = 1 for |m — n| = 1, Xmn = 1 
for |m — n| — V2 and Cnm — 1 for |m — n| = 2. The values A", B", C" and Z?" arc listed in appendix 0. The 
Fourier transformation is now a simple calculation. The approximation for small wave vectors (continuous limit) is 
by a special interest. The actual lattice structure becomes irrelevant on these sufficiently large spatial scales, i.e. the 
Fourier transformed are isotropic values. One obtains up to the second order in q: 

r"('^)-^(Z^(ro-7V) (32) 

The coefficients Fq and 7" follows immediately from the values A", B", C" and 13": 

Ff^ = A" + ZcS" + izc(^c - 2)C" + Zc-D" and 7" = + (z^ - 2)C" + 4L>" (33) 

(zc is the coordination number of the d-dimensional lattice, i.e. is the number of nearest neighbors per lattice cell. 
Straightforwardly, the Fourier transformed matrix g"^{q) can be written as 

with W = Tq — (FJ)^ — T^?^- Finally, the frequency matrices Q^^^^q) and Q^^^\q) are given by 

- ^ - + {2rh' - rj7^) q'} (35) 

and 



7vi7(")(9) = {rg - rlrl + (fJ7^ - 7^) q'} 

in the continuous limit. Note, that (p3), (pq) and (Bq) considers already that 7" = 7^ = and Fq = 1. 



(36) 



V. ANALYSIS OF THE RELAXATION BEHAVIOR 



A rough understanding of the results so far is possible by an analysis of the relaxation behavior of the correlation 
fimction $(q,<) neglecting the memory terms containing in (|27|). In this case the correlation function is reduced to a 
finite continued fraction. It can be expected that this case is related to the high temperature limit corresponding to 
a more or less exponential decay of the correlation function. Furthermore, this approximation should be reasonable 
for the description of the correlation function at short time scales. Note that because of if^^"-' (g, i) ^const. and 
K^^^^{q,t) — ^const. for t — > 0, one obtains K^^^\q,z) ^ z^^ and K^^^'>(q,z) ~ for z ^ 00. Therefore, the 
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memory terms can be neglected at sufBciently short time scales t ^ or z ^ oo. The initial conditions $0(9) and 
$0(9) of the correlation function $ arc defined by equilibrium averages: 



*rmi(0) = (flnflm) = ^nm and $,im(0) = (fl,nLi]n) = (37) 



and consequently 



.o(.) = i44 = ivr^(.) 



(38) 



The normalized correlation function $(q, z) = $(g, 2;)/$o('z) follows from ( |27D under consideration of (^8|) and under 
neglect of memory terms. A simple calculation leads to 

$(o z) = ^ + Nn('^H<l) + NTHq) 



which can be written as 



with the poles 



and the intensities 



^1/2 



$q,z = + ^— 40) 

Z — Zl Z — Z2 



Arf](ii)(g) =p ^ {NVti^^){q)Y - 4iV17(io)((?) . (41) 



jVri(g)-^2 , ^i-A^ri(g) 
= A2 — 

Zl - Z2 Zl - Z2 

The present approximation of the SFM[2, d\ is characterized by two relaxation times Tj^iq) = z^^ and T^{q) = z^^. It 
can be verified by a simple calculation that both relaxation times shows only a weak dependence on q. Furthermore, 
both relaxation times approach finite values for g — > 0. Obviously, the spin facilitated kinetic Ising model shows no 
diffusion-like modes which behave as r '--^ q~^ for the limit g — > 0. This finding agrees with investigations of the one 



dimensional spin facilitated kinetic Ising model |33|. The SFM[1, 1] corresponds to diffusion processes combined with 
creation and annihilation processes of active states. A cell has an active state if this cell can change its state without 
any support by further flip processes of neighbored cells. Creation and annihilation processes dominates at sufficiently 
large scales, i.e. an inhomogeneity reaching over a sufficiently large distance will be reduced by local creation and 
annihilation processes of mobile cells until diffusion processes becomes effective. 

As above mentioned, we restrict our investigations to the borderline case of macroscopic scales, i.e. g — > 0. Fig.|l| 
shows the relaxation times t^((z) and T^{q) as a function of temperature. As expected, there is no significant difference 
between the relaxation times for q = and g ^ 0, respectively. Furthermore, the slow relaxation time T^(g) shows a 
weak non-Arrhenius behavior. One obtains lnT|, In v{T) + o(q^) (see also eq.|^) and Inr^ '-^ In v{T') + 2£/T + 0(9^), 
respectively, at low temperatures. On the other hand, the high temperature regime is characterized by another 
temperature dependence: Inrj'j'^ = Ini^(r) + u^-^e/T + o(g^). The coefficients u^'^ depend on the actual lattice 
structure and the spatial dimension. But a simple analysis shows that always < 2, i.e. the activation energy of the 
SFM[2, d] increases with decreasing temperature. 

The existence of two relaxation times means not that the SFM[2,d] is characterized by an a- and a /3-process. 
The superposition of both decays, exp {— t/Tj^} and exp {—t/T^}, considering of the intensities Ai and A2 (see fig.|]), 
shows a continuous decay of the correlation function $(g, t), see fig.^j. This behavior is in an agreement with numerical 
simulations (l^Jl^ and it corresponds also to the above discussed thought, that spin facilitated kinetic Ising models 
are possible candidates modelling the behavior of supercooled liquids below the critical temperature of the usual mode 
coupling theory. 



VI. DETERMINATION OF THE MEMORY MATRICES 
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A. Complete and orthogonal basis 



All operators acting on the Fock-space can be represented by a complete collection of orthogonal base operators. 
The determination of such a basis is possible under consideration of the underlying di, (pseudo fermionic) algebra 
of the SFM[2, d] . The base operators can be expressed as all possible products of the above introduced operators % 
(^. A base operator is denoted as B^'' . (The index n corresponds to the order of the product, N„ is an n-dimensional 
vector indicating the concerning lattice cells). The first groups of the basis are: 

m = 1 

B^f = m-fij for i < j 
Bijk = ViVjVk for i < j < k 

Note that because of the commutation relation — the components of N„ can be ordered. The case of two 

or more equivalent indices is excluded because rjf = {1 — 2CToq)/CTeq(l — 'ffeq)f}i + 1, i-e. quadratic or higher powers of 
each operator f]i can be always reduced to a linear representation. The orthogonality means that: 

b(;)b^™))=5""Mn„,n„ (42) 

This relation can be checked considering that all equilibrium averages of operators on various cells decay in a product 

of averages with respect to these cells, e.g. {fj^...f}j fii----) — {vf) ■■■ iVj) {vl) This important relation is 

valid for all SFM[n, d] if the neighbor-neighbor interaction vanishes. Note that the Hamiltonian of the analyzed class 
of spin facilitated kinetic Ising models is given hy H = eDi{t) ~ - ey^^cxj(l^^cq)%(0- From this point of view, 
the relation ( ^2| ) follows immediately because of (?);) — and (t)!) = 1. Thus, the basis B = |-Bn^} orthogonal. 
The completeness of S = understood in relation to the reference state (s|, i.e. the following equation 

is fulfilled for an arbitrary operator X: 

(^l^ = EE(^^Ni)(^|S^"i (43) 

n N„ 

The mathematical proof of this property is given in appendix Bl 



B. Decomposition of the memory terms 
Eq.|l8| is a reasonable starting point for an analysis of the memory terms. The consideration of (^) and (^ leads to: 

^S"^ (^) = E {vkL-'^'Qexp {QLQt} QL'^^^^ gf (44) 

(Note that ~ Q). The apphcation of the completeness relation ( ^3| ) onto ( [4^ ) yields: 

4f W = EE E <Ni(<exp {QLQt} s(;:))H(rj -5^^^ (45) 

with 

H:il = {flkP'-'QB^^l) and H^!l = {B^^lQL-'+\lk) (46) 

These coefficients can be determined by simple algebraic calculations. One obtains immediately that both, -ff^^^^ ^^'^ 
H^^" vanish identically for n = 0, 1. On the other hand, the coefficients (46) vanish also identically for n > 5 because 

7 < 1. Hence, the memory terms can be constructed by using a finite number of functions ^^^^ ^-'^P |Q-^Q^| ^n'^^ 
These functions will be transformed identically. One obtains: 
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(b^"^ exp [QLQt] B^,";)) = (b^,") exp [QLQt] exp {-Li} exp {Lt} 

= Ep En, (^n^ exp (qLqA exp i-Lt\ B^^l") /b^^I exp {it} B^Tj) (47) 
= Ep En, exp {oLQt} exp {-it} B^,^;) (s^,^i(t)<i) 

The averages /B^\t)B^'^ \ are usual many point correlation functions. 



C. Mode coupling approximation 

1. Short time evolution of the memory 

An exact determination of (^B^^ exp {QiQt} exp {— ii} B^^J^ and ^i3^^(t)i?^"^^^ may be possible only for some 

few special cases. Therefore, we need a suitable approximation for a further treatment. In a first step we analyze the 
function: 

^nTn.W = {<lexp{QLQt}exp{-Lt] B^^l) (48) 

(^) can be interpreted as a separation of fast and slow time scales. The operator QLQ is related to a dynamics 
completely different to the dynamics of L. In general, it can be expected that 5^-* exp shows a significant 

evolution on a very short time scale in comparison to the characteristic time scale related to B^^ exp {^^}- Therefore, 

we come to the rough conclusion: while the evolution operator L contains all relevant time scales, the operator QLQ 
is mainly determined by contributions related to short time scales, i.e. the long time contributions are partially 
cancelled by the projection procedure. Consequently, the term exp {qI/Q*} exp {— it} and therefore ^'^''n (0 
should be dominated by long time scales, because the fast time scales are eliminated at least partially by the factor 
exp{— Lt}. Clearly, this is a very rough interpretation, but it gives an explanation for the assumption that the 

time dependence of ^^^^ (t) is weak in comparison to the decay of the correlation function (^B^\t) B^'^'^ which is 

connected only with the time evolution factor exp {^^}- 

However, it seems to be reasonable to expand ^'^''n (^) powers of the time t. The determination of all Taylor 
coefficients is out of the question. But an approximative analysis is possible by using a finite number of coefficients. 
We hope that the error of this perturbation expansion is sufficiently strong suppressed by the corresponding factor 



B^\t)B^j^J) (seeeqg7|). One obtains: 



°o .M 

*N„Npl''J - ^>-N„,Np ^1 

M=0 



with: 



. (np),M _ / A(n) d 
^^N„,N„ - \ ^N„ g^M 



Mr. , r . 

(p) 



exp {(5L(5t}exp{-Lt} 



Bjsj 



The first coefficients A^''"']^^ can be determined by simple calculations, e.g. 

.(np),0 _ / r(") r(p)\ — 
^'-N N — \^N ^N / ~" 0N„,N„ 

/ .X . \ (49) 

^Nr,Np = \-Sn„ [Q 'LQ - Lj B^^j 

In principle, the discussed expansion is an exact representation (The radius of convergence of the exponential function 
is infinite large). The true approximation consists in the breaking of the Taylor expansion after a finite power of t. 
We restrict our investigation to the simplest case, i.e. we assume A^'^)^"^ = for M > 1. Hence, one obtains 
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*i.TN,W«^"''^N„,N, (50) 

But it should be remarked that an extension to higher terms is possible without any problems. We abstain from a 
consideration of higher terms with respect to the clarity of the calculations. Furthermore, the obtained results (see 
below) using (|50|) show already a reasonable agreement with numerical simulations. 



2. Decomposition of the many point correlation functions 

The main problem is a reasonable approximation of the function {^B^"* (i)i3^*'' ^ This function decays in products 

of simple pair correlation functions if the distances between the corresponding lattice points (defined by the vectors 
Np and N^) are sufficiently large: 

- ^ ($,,,(t)$.,,,(t)....$.,,,(i) +perm) (51) 

This asymptotic limit is correct for infinitely large (or at least sufficiently large) distances between the lattice cells 
«i, 12, .... We use this borderline case as an approximation for an arbitrary set of lattice cells {Np,Nm}. This 
approximation is equivalent to the decomposition of higher static correlation functions into simple pair correlation 
functions. For example, a similar approach was used to create self consistent equations for the static structure factor 
p^ . Furthermore, this approximation is also the kernel of the well known mode coupling approach ||35[| . 



3. Reduction of the basis 

The third approximation consists in a reduction of the basis B ~ |^n^|- demonstrated that only base 

operators with n < 5 are necessary for a representation of (^6|). Because of the approximations (|5^) and (|l]) all 

other factors of (|4|) contains also no higher base operators . The consideration of all relevant is no general 
problem. However, we restrict the calculations only to elements with n < 2, also with respect to clarity. Really, the 

" (n) 

neglected terms containing base operators with 5 > n > 3 yield only small additional contributions to the final 
results. Thus, the complete representation of the memory ( p5[ ) is given by the approximation: 

a,i iii2jlj2 



D. Macroscopic scale 



As above mentioned, the g-dependence of the frequency matrices is very small. This weak dependence can be 
expected also for the memory terms, i.e. we restrict out investigations only to the macroscopic scale g — *■ 0. The 
assumption of a weak dependence on q leads to: 



^nm(i) = 'P{t)5n,n = ^ E '^^^^ ^ m)} 



(53) 



This approximation is related to the fact that all correlation functions between different lattice cells $,mi(i) with 
n ^ m (see eq.^) vanish for t = (<i>„m(0) is the equilibrium average {fjnfjm) which decouples as a result of the 
simple Hamiltonian H ~ '^if}i{i), i-e. {finflm} = (fin) (flm) = 0) and for t ^ oo (because of the ergodicity follows 
again $„m(oo) = (?)n(oo)r7m(0)) — (77„(oo)) (f?m(0)) = 0). Furthermore, one obtains only a very small correlation 
^nm(t) between different cells for finite times t differences which can be checked by numerical investigation. 
From (|38| ) and ( ^2| ) it follows for g ^ the relation (7o(0) = — FJ/tq. Thus, the evolution equation on a macroscopic 
scale is given by (Note, that there is the initial condition: (p{0) = 1): 



ip{z) 



iVr!(io) (0) - iVi^(io) (0, z) + zTI/to 



z + 7Vf^(") (0) - Aris:(ii)(0, z) - ri/ro 



(54) 
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The determination of the frequency matrices was reaHzed above. One obtains in the macroscopic hmit by using (|3£ 
and M): 



^^''''W-y-prM- and iVl^("HO) = f^F^ (55) 



It remains the determination of the memory. Using the definition: 



and (63), it follows from (B 



h"ici,ci') = ^E,,M5S)exp{*q(k-j) + *q'(l-j)} 
/i"(q,q') - ^E,,M^2)" exp{*q(k-j) + iq'(l-j)} 



a 

« X;^''(0,0)¥'(t)2/i"(0,0)7V.9"^(0) 

a 

The substitution ^.'''(q, — q) h^{0,Q) and /i"(— q, q) — > ft-"(0,0) is possible because the g-dependence of these 
quantities is again relatively weak. We need /i^(0, 0), /i''(0, 0) and h^{0,Q) for the following investigations. A simple 
calculation leads to /i°(0, 0) = 0. Using (^) the memory terms can be written as: 

iVi^(io)(0,t) = Toip{t)^h\0,0)h\0,0)^^^^ 
JVX(")(0,t) = ToV(^)'feHO,0)^HO,0) rg-(r'ip 
It should be remarked that the ratio between both memory terms is given by the relation: 

(0,i) ^ TO 

Finally, we must determine the quantity A = Tq /i^(0, 0)/i^(0, 0)/(4(rQ — (Fq)^)). This is again an algebraic procedure 
which can be realized straightforwardly. The final results are very unwieldy. Therefore, we give only the explicit 
expressions for the asymptotic case of low temperatures, i.e. for acq 0. One obtains for a square lattice (zc ~ 4) 
and a cubic lattice (z^ = 6), respectively: 

A = ^(l + o(aeq)) for Ze = 4 
A = ^(l + o(aeq)) for = 6 

The behavior A ^ 5"^^ is characteristic for T — s- and a^cq — 0, respectively. (^) can now be written as: 



Tl 1 Fi7Vf](ii) (0)ro - iV17(i") (0)r2 - (F^) ' 



TO TO zTo + iVri(") (0)to - XE{z) - Fi 
with S(z) = (dt/T-o) (^(0' exp{-zO. 



(56) 



VII. DISCUSSION 



Now we are able to analyze the characteristic slowing down of the dynamics of the SFM[2, d] for decreasing tem- 
perature. The central equation for the following discussion is (|56|). The first question is the existence of ergodicity 
and nonergodicity. Exists there a critical temperature T*, so that the correlation function ip{t) shows an incomplete 
decay ip{t oo) = foo ^ for T < T*? In other words, has the function (p{z) a pole at z = below the critical 
temperature T*7 This question is equivalent to the determination of a kinetic phase transition from an ergodic state 
into a nonergodic state for supercooled liquids . To this aim we split the correlation function into a nonergodicity 
part and a contribution ipcrg{t): 
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<f{t) = /oo + </5crg(<) 



(57) 



The function (^crg(0 describes the remaining ergodic part of the SFM[2, d], i.e. fcig{t oo) = 0. The Laplace 
transformation leads to 

viz) ^ ^ + -^ci-giz) (58) 
with lini zip crg{z) = 0. The memory term S(z) can expressed by: 

S(z) = ^ + Serg(z) (59) 

ZTq 



with lini zSerg(2) = 0. 1^ , (|59|) and (^ yield 



lim zip{z) = foo = lim 

z—*0 z^O 



^ ^ Th ^ riiV»(")(0)ro-iV»(i°)(0)ro^-(ri)^ 

ZTQ A/^ 



(60) 



It follows immediately that the nonergodicity part /oo has a nonvanishing value only if = . Otherwise, the only 
solution of (|6^) is given by foo = 0, i.e. the SFM[2, d\ is an ergodic system if Fq ^ 0. The value of Fj vanishes only 
for T = 0: 

Fj = T = 

i.e. the SFM[2, d] realizes a kinetic phase transition from an ergodic system to a nonergodic system at the critical 
temperature T* = 0. Thus, the only nonergodic state can be observed at zero-point temperature. Additionally, we 
obtain from (M): 



rj(")(0) = and f7(i°)(0)=0 for T = 
Thus, the nonergodic part is given by: 

/oo = 1 

i.e. an initial equilibrium configuration at T — shows no structure relaxations during the total observation time. 
This behavior is a consequence of the fact that the kinetic phase transition occurs at the absolute zero at temperature. 
In other words, each arbitrary equilibrium configuration is frozen at T = = 0. 

We obtain the important result that the SFM[2, d] is ergodic for all finite temperatures T > 0. This is a real 
contradiction to the statements of the original papers p^ , p^ which predict a kinetic phase transition at a finite 
critical temperature. 

Of course, it is possible to explain T* = by a simple picture. A finite fraction of active cells (i.e. cells which are able 
to change their state without any previous change of states of neighbored cells) exists at each temperature T > 0. 
The concentration of these active or nonfrozen cells is proportional to a^^ (a^q is equivalent to the probability that 
a cell has two neighbored cells with cr = 1). Active cells are mainly isolated at sufficiently low temperatures, but 
an annihilation of an isolated active cell is never possible. The property to be an active cell can be transferred to 
a neighbored cell by some few elementary flip processes, (this procedure can be interpreted as a diffusion of active 
states), a new active cell can be created in the nearest environment of an initially active cell (creation process) and 
two neighbored active cells can be unified to one active region and further this region can be reduced to one active cell 
(annihilation process). Diffusion, annihilation and creation are well-balanced in the thermodynamical equilibrium. 
On the other hand, this three processes realize a motion of active states through the whole volume, i.e. each cell is 
able to change their state after a sufficiently long waiting time. Note that this considerations must be modified for a 
finite volume because some configurations of the SFM[2, d] show a self blockade in finite geometries |1^. But it should 
be denoted also, that a special choice of the initial configuration excludes any type of self blockades |l^]. However, 
the main result of our discussion is the ergodicity of the SFM[2, d] for T > 0. 
To analyze the relaxation behavior at finite temperatures T > we introduce the relaxation time 

_ _ iV»(")(0)ro-AS(0)-Fi 

^° iVl](io)(0)r2- AFiS(O) ^ ^ 
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and the coefRcient 



^ = Arr!(")(0)To-AS(0)-rJ 



(62) 



Using these notations, we get from (pq): 



ip{z) = 



To 



zTa + A [S(0) - 



(63) 



(IsTI), ( |62|) and ( |63| ) are a closed, nonlinear system of equations, which can be solved by numerical standard methods. 
Fig.^ shows ip{t) for a SFM[2, 3] on a square lattice {zc = 6) for various relative temperatures T/e. We see that the 
correlation function Lp{t) shows with decreasing temperature a pronounced stretched decay over some decades, while 
an exponential-like decay is obtained for high temperatures. This stretching can be illustrated by a simple argument. 
Short times (t ^ or 2; — > 00) are related to a behavior Lp{z) ~ (z + Fj/ro)^^ or ip{t) ~ exp {— Fjt/ro}. On the 
other hand, the long time regime (t ^ 00 or z ^ 0) is characterized by ip{z) ~ (z + t~^)~^ or ip{t) ~ exp{— t/rc}. 
As there is FJ/tq » we expect a typical crossover between both regimes characterized by a stretched decay, see 
figj^. But it should be remarked, that the decay of ip{t) is no Kohlrausch-WiUiams-Watts function (^ exp(— (t/r)^). 
There exists only over a finite interval a reasonable fit with such a stretched exponential function [ p7[ . 

The spectral density S{Xs) of the correlation function is defined as the set of amplitudes of exponential decays 
which contribute to (p{t). Thus, ip{t) is the Laplace transformed spectral density with respect to the Laplace variable 
t: 



Lp{t) = / dA,S'(A,)exp(-A,0 



(64) 



One obtains the remarkable result that the spectral density is positive definite, (fig-^). The knowledge of the spectral 
densities allows the determination of other interesting properties, for example the susceptibility x{'^)^ see fig. |6|. 
Finally, the averaged relaxation time t(T) can be obtained by using 



-(T) 



/ tip{t)dt J S{Xs)X-'^dXs 




/ ^it)dt J SiX,)Xi^dXs 




(65) 



This relaxation time shows with decreasing temperature an increasing deviation from the simple relaxation times 
(|l|). Especially, one obtains now a typical non-Arrhenius behavior (fig.|^) . 



VIII. CONCLUSIONS 



It was shown that irreversible master equations can be easily transformed into projection equations by using the 
Fock space representation. Whereas the usual projection formalisms, which leads to the well known Mori-Zwanzig 
equations, start from a reversible Liouville equation, the master equations are already irreversible. As a consequence 
of this initial irreversibility, one obtains additional frequency matrices and memory terms. Thus, these additional 
contributions are caused mainly by the loss of the invariance against an inversion of the time. 

A second important property follows from a general analysis of the frequency matrices. The corresponding poles of the 
Laplace transformed correlation function $((7,2;) are always located on the negative real axis. Especially, the imagi- 
nary part of the poles vanishes identically. This behavior is related to the general structure of the master equation. 
The dynamical matrix L{a, a') of a master equation is always negative definite (or better semidefinite because at least 
one eigen value is zero as result of the conservation of the probability). Thus, only relaxation processes should be 
observed, i.e. the evolution of the probability P {a, t) can be approached by a probably infinitely large expansion in 
terms of exponential functions: P {a, t) = P^q {cr, t)+J2m ^rn (c) exp {—k^t] with (A^ > for all m). One obtains no 
oscillations in contradiction to microscopical systems which bases on a Liouville equation. The absent of oscillations 
is directly connected with the Markov property of the underlying master equation. This behavior corresponds also to 
the fact that the observation of a spin wave propagation (corresponding to density waves in a glass or a supercooled 
liquid) is not possible for the SFM[n, d]. From this point of view, the traditional notation 'frequency matrix' can 
be misleading. It seems to be possibly favorable to use the notion relaxation matrix. However, we have used the 



14 



traditional terminology to avoid misunderstandings and conflicts with other well defined quantities. 
In principle, the consideration of all possible time derivations fif\t) = d^f]i(t)/dt^ = f]i(t)L^ (/3 = 0,1, ...,oo) as 
relevant operators leads to an infinite continuous fraction for the correlation function $„m {t) , determined by frequency 
matrices of various order. This representation allows a systematic analysis of the short time behavior of the SFM[2, d], 
because an infinite continuous fraction contains no memory term. Unfortunately, a general explicit determination 
of the frequency matrices, e.g. by a successive rule, cannot be obtained. On the other hand, each approximation 
using a finite number of frequency matrices (e.g. fi^^"' and fi^^'^' in the present case) fails for sufficiently long times. 
However, some important properties of the spin facilitated kinetic Ising model can be verified for this relatively rough 
approximation, e.g. a weak non-Arrhenius behavior of the relaxation times. Note, that this result is in agreement 



with various numerical simulations [g^ 19 1. On the other hand, the typical stretched decay of the correlation function 
can not be explained by using this simple approximation. 

A satisfactory treatment is possible by an approximative consideration of the memory terms. This procedure leads to 
an equation, which is partially similar to the well known mode coupling equation of supercooled liquids. The memory 
terms yield the main contributions to the typical stretching of the correlation function. Some small remaining devia- 
tions from the numerical results are caused by the approximations ( ^0|) and ( ^l|) of the memory terms. 
From this point of view, we come back to the initial question. Which processes of the usual glass transition can be de- 
scribed by the spin facilitated kinetic Ising models? Obviously, the spin facilitated kinetic Ising model is not adequate 
for a description of the fast processes inside a supercooled liquid. This statement is supported by both, numerical 
simulations ||l9[| and the presented analytical investigations, which show that no fast processes can be observed. 
On the other hand, spin facilitated kinetic Ising models allow a more or less reasonable, quantitative description of 
the slow (a-) process in supercooled liquids below the critical temperature of the usual mode coupling theory [|l|,Q . 
Note that for this low temperature regime the time scales of a- and /9-process are well separated, i.e. a separation 
of the underlying dynamic is actually possible and it is considered in the structure of the master equations related to 
the SFM[n, d], see above. The fast dynamics, corresponding to the /3-process, determines the thermodynamical noise. 
This noise is not explicitly contained in the previous equations, but it is the underlying cause for the irreversibility 
of the master equations. Additionally, master equations and stochastic evolution equations are principally equivalent 
p6[ . However, the remaining slow dynamics of a supercooled liquid (a-process) is represented in the kinetic scenario 
of the SFM[n,d]. 
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APPENDIX A: REPRESENTATION OF T 



The rigorous values A", _B", C" and D" of ( ^3| ) can be obtained straightforwardly by using the evolution operator 
L the commutation relations (^) and the definition of F^'j, (29). Under consideration of the coordination number 
Zc and the thermodynamical equilibrium acq of the cell state we obtain the following results ^"-terms: 



_B"-terms: 



^0 

A^ 
A^ 



-2 

oq 

2 



4a: 
8a 



Zc 

2 

2 

Zc 

2 



6o-o 



24cto 



6< 



180ct^ 



+ Aai^ (1 - CTcq) 2;c(Zc - 1)^ (1 + CTcq(2c - 2)) (2 + CTcq(Zc - 4)) 



90<y 6 



BO 



4^3 (l-aeq)(Zc-l)2 

5=* = 16a^q(l-aeq)zc(^c-l)'(l + 3^cq(2c-2) + CT2q(z,-2)(z,-3)) 



C" -terms: 



and 



= = = 



- 8a3q(l - aeq)2 (2 + 2acqZc{Zc - 2)) 



Z3"-terms: 



and 



= = Z?^ = 



^a^q(l-aeq)2(^e-l)' 



APPENDIX B: COMPLETENESS OF THE BASIS B 

The proof consists in two parts. At first we analyze an operator X' , which is a multilinear form of the operators fji 
i.e. 



(Bl) 



N„ 



with arbitrary coefficients /^n"*- The orthogonality of the basis B = l-B^'j leads immediately to /J^^' = l^X'B 
i.e. all coefficients of X' can be determined by a successive procedure. In other words, the basis I-Bn"*} complete 



with respect all operators X', i.e. the basis forms a space containing all operators of type (Bl). 



The second part of the proof analyses arbitrary operators X acting on the Fock-space. X consists in operators f]k as 



well as annihilation operators dk and creation operators c?J,. A representation like (Bl), extended by the operators d, 
and d\ is always possible: 



X 



Mi.M2,M3.M'' k=l 



N p 



(B2) 
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(i is the simple unit operator 1 = 1, A'^ is the number of lattice cells). The vectors M''' (7 = 1, ...4) contains N integer 
numbers, i.e. = {mj, ...,mj^} with > for all possible k and 7. The commutation relations (||) can be used 
to write X as a representation with the internal restrictions: mJ, + + m| + mj. — 1 for all lattice cells k, i.e. each 
contribution to the sum (B2) contains exactly one of the four operators i, dk, dj, and fjk with respect to any cell k. 
Furthermore, the commutation relation (||) allows a shift of all operators related to a given lattice cell i to the left 
hand side. It is simple to show that P8|: 



= (s\d,^(s\ l-a. 



(s| % = (s| % (s| = (s| [acq + fliy^acqil " CToq)) 



Obviously, the application of the operator X onto the reference state (s| is equivalent to the application of a corre- 
sponding reduced operator X' containing only operators % (and the trivial operators 1 = 1) onto the state (s|, i.e. 
there is a definitely mapping 



X ^ X' 



with 



(six ^ (six' 



Therefore, one obtains (see eq.||): 



F 



X'B^^l 



F 



X'B^^l 



and consequently by using (Bl) 



.1 x^(s\x'^Y.ll ) (^1 < -EE ) (^1 B. 



n N„ 



n N„ 



i.e. an arbitrary operator X of the Fock-space can be completely presented by using the basis B — |-Bn'} 
the consideration that this operator acts into the left direction on the reference state (s\. 
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FIG. 1. Arrhenius-plot of the relaxation times r^(g) and r^(g) as a function of the reduced inverse temperature e/T. The 
full lines correspond to q — 0, the dashed lines are related to q = {I is the unit size of the lattice cells). The dotted line 
corresponds to an Arrhenius like process. The deviation between the dotted line and the full line indicates a weak non-Arrhenius 
behavior of T^(g). 

FIG. 2. Intensities Ai and A2 as a function of the reduced inverse temperature e/T. The full lines correspond to g = 0, the 
dashed lines are related to q = . 

FIG. 3. Correlation function "l>(g, t) for various reduced temperatures T/e = 6.0, 4.0, 2.0, 1.0, 0.5 and 0.1. The temperature 
decreases in the direction of the arrow. 

FIG. 4. Decay of the reduced correlation function (fi{t) for various reduced temperatures T/e — 0.11, 0.15, 0.22 and 0.45. 
The temperature decreases in the direction of the arrow. The functions show with decreasing temperatures an increase of the 
stretching. 



FIG. 5. Spectral density for various T/e = 0.11, 0.15, 0.22 and 0.45. 

FIG. 6. Suzeptibility x('^) for various reduced temperatures T/e = 0.11, 0.15, 0.22 and 0.45. The temperature decreases in 
the direction of the arrow. 



FIG. 7. Relaxation time t as a function of the inverse reduced temperatures e/T. 
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